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Compact directed percolation is known to appear at the endpoint of the directed percolation crit- 
ical line of the Domany-Kinzel cellular automaton in f + f dimension. Equivalently, such transition 
[•"t . occurs at zero temperature in a magnetic field H, upon changing the sign of H, in the one-dimensional 

Glauber-Ising model, with well-known exponents characterising spin-cluster growth. We have in- 
' vestigated here numerically these exponents in the non-equilibrium generalization (NEKIM) of the 

Glauber model in the vicinity of the parity-conserving phase transition point of the kinks. Criti- 
cal fluctuations on the level of kinks are found to affect drastically the characteristic exponents of 
spreading of spins while the hyperscaling relation holds in its form appropriate for compact clusters. 
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In the one-dimensional Domany-Kinzel automaton (DKCA) jl]|^] the state o-(i,t) of site i at time t depends on 
a(i — l,t — 1) + a(i + 1, t — 1), (er(i, t) — 0, 1). Of the conditional probabilities p{cr{i — 1, t — l)er(i + l,t — l)\a(i, t)) 
the independent ones are denoted by po = p(00|l), p\ = p(01|l) = p(10|l) and P2 = p(ll|l). All sites are updated 
simultaneously in the process. The phase diagram of the DKCA in the (pi,P2) plane, exhibits a line of (second order) 
critical points of directed percolation universality class, which line ends at the so called compact directed percolation 
point (CDP). This point is situated on the line pi = l,po = at pi = 1/2. By crossing this point (changing the sign 
of pi — 1/2) the transition is a first order one between two ordered phases (empty and full or, equivalently, using the 
spin variable s(i, t) — 2a(i,t) — 1, all spins up and all spins down). The characteristic critical exponents of the CDP 
transition are known exactly and a hyperscaling relation has also been derived for such transitions in d dimensions 
£^ ■ JUI)- For t ne spreading process of a single a(i, 0) = 1 in the sea of zeros the exponents S s , rj s and z s defined at the 
' transition point for the power-law time-dependences of the density of l's ti s cx f 1 ' , the survival probability P s (t) oc t~ Ss 
t— I ■ and the mean square distance of spreading < R 2 s (t) >~ t Zs have been obtained as 0, 1/2 and 1, respectively ||. (In 
the following subscript s refers to spins for all the quantities. Without subscript the corresponding quantity for kinks 
is meant, except for v). For the parallel(time-direction) and perpendicular (space-direction) coherence lengths and 
v±_ resp., as well as for the dynamical critical exponent Z Domany and Kinzel has obtained the exact results: v\\ = 2, 
v± = 1, Z = 2, respectively (which means only two exponents as by definition v\\ = Zu±_). The above mentioned 
hyperscaling law H 



T] s +S s = dz s /2 (1) 



is fulfilled with the above exponents. More generally, Dickman and Tretyakov argue, that eq.([l]) is valid at first order 
£j , transitions and it should apply to cases, whenever power-law growth produces compact 'colonies', developing from 
single seeds. In ref. Q 'compactness' is clearly defined: it is meant that the density of colonies in surviving samples 
. ' remains finite for t — > oo. 

It is obvious, that the above sketched 1+1 dimensional CDP transition is equivalent to that in the 1-d (ferromagnetic) 
Glauber-Ising model Q at T — 0,because the symmetry as well as the kinetics are the same. Changing the parameter 
Pi of the DKCA around p\ = 1/2 corresponds to introducing a magnetic field H into the spin-flip probability wi of 
the Glauber-Ising model (Section II.) and changing its sign. 

On the basis of this equivalence it is of some interest to investigate the same spreading problem in the framework 
of the nonequilibrium generalization |^,^] of the kinetic Ising model (NEKIM) (Section III.), where, in some range 
of its parameters, there is a continuous transition between a single-domain and a multidomain state. The order 
parameter of this transition is the density of kinks. The critical fluctuations of this so called parity conserving (PC) 
transition 0- nave pronounced effect on the underlying spin system, as was found earlier by the authors 

p2| both in case of static and dynamic exponents in situations of quenching from T = oo (random initial states). 
These investigations will now be completed by studying, via numerical simulations, the spin spreading process at the 
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PC point (Section IV). It is found that the characteristic exponents differ from those of the CDP transition, as could 
be expected, but still basic similarities remain. Thus the transition which takes place upon changing the sign of the 
magnetic field is of first order and its exponents satisfy eq. ([!]). Accordingly it can be termed as 'compact', and we 
will call it compact parity conserving transition (CPC).Q The static magnetic critical exponent A is also determined 
at the PC point. 

II. GLAUBER-ISING MODEL 

The d = 1 Ising model with Glauber kinetics is exactly solvable. In this case the critical temperature is at T = 0, 
the transition is of first order. We recall that px — e^^r plays the role of T ^ r Ta in Id and in the vicinity of T = 
critical exponents can be defined as powers of px, thus e.g. that of the coherence length, ^, via f oc px" ■ In the 
presence of a magnetic field if,( when the Ising Hamiltonian is given by Ti = —J^2 i SiSi + \ — HY], Sj, s, = ±1), the 
magnetization is known exactly. At T = 

m(T = 0,ff) = sgn(H). (2) 
Moreover, for £ 3> 1 and H/kT <C 1 the the exact solution reduces to 

m~2ft£; h — H/ksT. (3) 

In scaling form one writes: 

m-C^g(h^) (4) 

where A is the static magnetic critical exponent. Comparison of eqs. (|J) and (|J) results in j3 s = and A = v . These 
values are the ones well known for the 1-d Ising model. It is clear that the transition is discontinuous at if = also 
when changing if from positive to negative values, see eq.(|2])( In the following the order of limits will always be meant 
as: 1). if -> and then 2). T -»■ 0). 

The kinetics of the Ising model in a magnetic field has been formulated by Glauber Q|. In its most general form 
the spin-flip transition rate for spin sitting at site i is : 

Wi h = Wi(l — tanh hsi) w Wi(l — hsi) (5) 

Wi = -(1 + Ssi-iSi+i) (1 - -Sifa-i + s l+1 )j (6) 

where 7 = tanh 2 J/kT (J denoting the coupling constant in the Ising Hamiltonian), T and 5 are further parameters. 
This model will reach the same equilibrium state as the Ising model in a magnetic field. 

For the case 6 = 0, T = 1.0, which is usually referred to as the Glauber-Ising model, Z = 2 (Z is the usual dynamic 
critical exponent) is also a well-known result. The non-equilibrium generalization of the kinetic Ising model to be 
used later on will be given a short review in the next section. 

III. THE NEKIM MODEL 

In the NEKIM model, besides the spin-flip transition-rate eq.Q, taken at T = 0, also a nearest neighbour mixing 
of spins with probability p ex is applied at each time step of the simulation. The spin-exchange transition rate of 



nearest neighbour spins (the Kawasaki |14| rate at T — 00) is lUij+i = |p ea; [l — SjSj_|_i], where p ex is the probability 
of spin exchange. Spin-flip and spin-exchange are then applied alternately. The model was originally proposed and 
investigated for values S > at finite temperatures in B. It is , however, at T = and for negative values of S, 
that in this system a second order phase transition takes place fl for the kinks from an absorbing to an active state, 



'in a previous paper of the present authors devoted to damage spreading investigations of different non-equilibrium 

one- dimensional models, the issue of a CPC transition has already been raised. 
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which belongs to the parity conserving (PC) universality class . The order parameter is the density of kinks, at the 
PC point it decays in time as a power law rikink oc t~ a , with a = .285(3). 

The absorbing phase is double degenerate, an initial state decays algebraically to the stationary state, which is one 
of the absorbing ones (all spins up or all spins down, provided the initial state has an even number of kinks) and the 
whole absorbing phase behaves like a critical point with power law decay of correlations, like the Glauber-Ising point 
(5 = 0,p ex = 0). 

Now let us look at the PC transition from the point of view of the underlying spin system. The above mentioned 
first order transition at T — of the Ising system disappears at the PC point and is, of course, absent in the whole 
active phase of the kinks. The fluctuations of this PC transition exert a pronounced effect on the underlying spin 
system as found earlier Jjjj thus e.g. the the classical dynamical exponent Z, defined, as usual through the relaxation 
time t s of the magnetization, t s oc £ z , was found to be Z = 1.75(1) instead of the Glauber-Ising value of Z — 2. In this 
case one approaches the PC point from the temperature 'direction', by decreasing it to (the effect of temperature is 
to create kink pairs inside of ordered spin-domains). On the other hand we can also look at the transition by changing 
a characteristic parameter (chosen by us to be 5) of NEKIM through the critical point S c and fixing the other two. 
As a function of e = (5 — 5 C ) | the transition - on the level of spins - is again a first order one of type order-disorder. 
Namely, taking initial states with an even number if kinks, the magnetization of the stationary state has a jump at 
e = 0. The same is true when changing the magnetic field h from negative to positive values at e = 0. Thus for the 
spins the value of the (static) critical exponent f3 s is zero, in all the three 'directions' of departing from PC (pt, e and 
h), as mentioned above. (For simulational results see |l2| , |l3| ]). 

In the following we will choose the same PC transition point as in previous works pSlJl^l, and make simulations at 
and around this point by changing the magnetic field h. The parameters chosen are: T = .35,p ex — .3, 8 C = —.395(2). 
In these previous simulations the spin-flip part has been applied using two-sublattice updating. After that we have 
stored the states of the spins and made L (L is the size of the system) random attempts of exchange using always the 
stored situation for the states of the spins before updating. All these together has been counted as one time-step of 
updating. ( Usual Monte Carlo update in this last step enhances the effect of p ex and leads to S c = —.362(1)). 



IV. SPIN-CLUSTER-GROWTH SIMULATIONS 



Time - dependent simulations have proven to be a very efficient method for determining critical exponents (besides 
the critical point itself) |l5|-|l7j]. On the basis of eq. ([!]), the f-dependence of the magnetization in scaling form can be 
written as 

m(t,h) ~ t~^g(ht^) (7) 

Such form can be used in a quench from T = oo to T c and was exploited also in ]l2] ] , though at h = 0, using temperature 
as a second variable, for determining mainly static critical exponents of the spins at the PC point. 

In the following we will further study the influence of the PC transition on the spin system from a different point 
of view. Instead of starting with an initial state of randomly distributed up- and down-spins with zero average 
magnetization as in the above mentioned simulations of quenching, we will now investigate the evolution of the 
nonequilibrium system from an almost perfectly magnetized initial state (or rather an ensemble of such states) . This 
state is prepared in such a way that a sin gle up-spin is placed in the sea of down-spins at L/2. Using the language of 
kinks (or particles, in the BARW model flqJl), this corresponds to the usual initial state of two nearest neighbour 
kinks placed at the origin. The quantities usually measured of the forming clusters are the order-parameter density, 
the survival probability P{t) and the average mean square size of spreading < R 2 (t) > from the center of the lattice. 
At the critical point these quantities exhibit power law behaviour in the limit of long times; more generally we can 
write 

n a (t,h)~fi-g l (hf&) (8) 
for the deviation of the spin density from its initial value, n s = m(t, h) — m(0), 

P s (t,h)~t- S *g 2 (ht&) (9) 

for the survival probability and 

<R s 2 (t,h)>~t z °g 3 (ht^) (10) 
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for the average mean square distance of spreading from the origin. The argument of the scaling functions above has 
been taken from eq.(0); but now at the PC point instead of the Glauber-Ising one. Thus the exponents A, v and Z 
in the above equations take values appropriate at the PC point. We note here that the coherence length exponent 
v appearing above is basically different from the v\_ and v\\ generally used in the context of DP transitions or in 
connection with the kinks in NEKIM. Namely, £j_ oc eT Vl - with e denoting the deviation from the PC point in the 
'direction' of the quantity driving the phase transition ((e =| 8 — 6 C | here). Moreover, i/m = i/±Z. Z is independent of 
the abovementioned 'directions', characterizes solely the transition point, as it should be jl2]. 

We have measured n s (t) and P s (t) at and in the vicinity of the critical point with initial configuration of a single 
up-spin at the origin in the see of down-spins and allowing the system to evolve according to the rule of NEKIM 
as described above. Averaging has been taken over runs with different sequences of random numbers during the 
evolution. 




1/t 

FIG. 1. Local slopes of the spin density n a (t) for zero magnetic field near the PC point. —5 — 0.393,0.394,0.395 (from 
bottom to top). The best scaling result is 773 = .288(4). In the averaging the number of independent runs was 3 — 5 x 10 6 . 
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FIG. 2. Local slopes of the survival probability P 3 {t) for zero magnetic field in the vicinity of the PC point. 
— S = 0.393,0.394,0.395 (from bottom to top). The best scaling result is S s = .287(3). In the averaging the number of 
independent runs was the same as for Fig.[j] 
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Fig.|l| and Fig. || show the local slopes (for definition see e.g. &) rj s and —6 S , respectively. As the survival probability 
must be the same for spins and kinks (if the minority spin dies out, kinks also disappear and vice versa) S s = S. The 
same applies also for the RMS size of the cluster. As no result has been reported before for 8 in the NEKIM model, 
exhibiting Fig.H has its own merits. r] s , however, is an independent new exponent. 
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FIG. 3. Phase diagram of NEKIM in the (h, — <5) plane. The chosen PC-point is at 5 — —.395. For lower values of —S, in the 
Glauber-Ising regime, the vertical line connecting the PC and the NGI points (at h = .0) consists of all CDP points with its 
characteristic critical exponents. The simulations around the PC point have been done here for h > in the interval < h < .1. 
The other parameters of NEKIM in the whole plane are as follows: V = .35 and p M = .3. 

(Besides results at the PC point we have also carried out detailed simulational studies at a point in the so called 
Ising phase, namely for 5 = 0, T = .35 and p ex = .3. This point is a non-equilibrium one due to the non-zero value 
of p ex , and is marked on Fig]3 with NGI(non-equilibrium Glauber-Ising) on the abscissa. The results which we have 
obtained via simulations at this point (Figj^) are, within error, the same as for the (exactly solved) Glauber-Ising 
case.) 

Fig.|| shows the the asymptotic values for large times of P s (t, h), for different values of h in the range of h = .005 — .1. 
For the exponent (3 J defined through 

lim P s (t,h) oc h p °' (11) 

t — >oo 

the value /V = .445(5) has been obtained. Fig.^ and FigJ^ show the scaling functions, eqs. (||) and @, respectively. 
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FIG. 4. The scaling function n s (t) oc f" and P s (t) <x t 
the data shown gives T] a NGI = .0006 and 8 a NGI = .500(5). 



3 at the NGI point of Fig.g ( p ex = .3, 5 = 0). The power law fit of 
Number of independent runs in the averaging was: (10 — 60) x 10 4 




FIG. 5. Level-off values of the survival probability P s (t,h) for large times at different values of h. The straight line (a) is 
power law fit near the PC point with f3 a ' = .445(5). Points around the straight line (b) show simulational results of the same 
quantity for the NGI point, giving f3 s ' — .99(2). Number of independent runs in the averaging was: 10 4 - 10 s . 
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FIG. 6. The scaling function n s (t,h). The different curves correspond to the following values of the parameter h: 
h — .05, .01, .009, .003, .001, .0005. The values of the parameters rj a , A and vZ are given in the text. Number of indepen- 
dent runs in the averaging was: 4 x 10 4 . 

Here the best fit for the scaling-together of data with different values of h could be achieved with A = .49(1), using 
the measured values S s — .285, rj s = .285, and that of vZ from former studies, vZ = .777 [H. 




FIG. 7. The scaling function P a (t,h). The values of the parameters rj 3 , A and vZ used are given in the text. Number of 
independent runs and values of h are the same as for FigJtl 
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Data for different values of h scale together sufficiently well when considering the relatively poor statistics ( averages 
over 4 x f 4 samples, typically). The scaling law 



/?.' = 



(12) 



following from eq.(^|) is satisfied with the above values of the exponents, within error. We note here that S ' can 
be connected with flkink = Pkink' , using in eq(|l2|) 6 S = 5 and the definition of fikink' from ref. | fl9f with the result: 

/V = f3 k inkv/(v±A). 

The hyperscaling law for the spreading exponents was derived in its most general form by Mendes et al. |pj| which 
we write here for the spin-quantities: 



(l + ^)5 s +ri s = dz s /2 

Ps 



(13) 



In eq. (|l3J) , in analogy with the spin cluster growth description at and in the vicinity of the CDP transition of the 
DKCA p[, the above finite value of S ' enters. (As explained in the introduction, (pi — p\ c ) of the DKCA, with 
Pi c — 1/2, corresponds to h in the Glauber-Ising formulation). Moreover, (3 S = 0, which value follows near the PC 
point from the same symmetry consideration as at the Glauber-Ising point (though does not in the active phase). 
Here, again, one should recall the above mentioned analogy between the DKCA's variable (pi — 1/2) and the variable 
h in the present case. With the exponents obtained and summarised on Table 1 eq.([l3j) is fulfilled. As already 
mentioned in the Introduction, according to the argumentation of || the fulfillment of the hyperscaling law in the 
above form is equivalent to compactness of the clusters. For illustration developing clusters are exhibited on Fig.^ 
under three conditions: A) Glauber case (CDP in the DKCA sense), B) at the NGI point (see Fig.||) where the kinetics 
is a non-equilibrium one (p ex ^ 0) and C) at the PC point. It is apparent that the minority phase never develops 
inside of the majority one, moreover, the branching process present in the kinetics in cases B) and C) makes the fla 
pieces of CDP boundaries fringed. 
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FIG. 8. Clusters developing from a single up-spin(dark point) in the sea of down-spins(white points) at t = for three choices 
of NEKIM parameters, see text 
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The results together with some of the critical exponents obtained earlier in ]12| are summarized in Table I. 
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.99(2) 


1/2 


.0006(4) 


.500(5) 


1(= 2/Z) 


CPC 


.00(2) 


.444 


.45(1) 


.49(1) 


.288(4) 


.287(3) 


1.14(= 2/Z) 



TABLE I. Spin-cluster critical exponents for NEKIM in a magnetic field 



V. SUMMARY 



In summary, we have carried out numerical studies of the power law behaviour of spreading of spins, at the PC 
transition point of NEKIM (where a second order transition occurs on the level of kinks). It has been found that 
the analogue of the Domany-Kinzcl CDP transition - a first order transition upon changing the sign of an applied 
magnetic field - still exists. Of the three exponents measured only A, which is the static magnetic exponent of the 
Ising model,was found to be unaffected by the critical fluctuations of the kinks within error. Concerning other static 
Ising exponents, this circumstance is not so natural, as e.g. the coherence length exponent is a counterexample (see 
ref. |T^] and Table 1). Consequently the relation A = v valid in the Glauber-Ising case is no more fulfilled at the PC 
point. f3 s ' characterizing the level-off values of the survival probability of the spin clusters is a new (static) exponent; 
A and j3 s ' are connected by a scaling law. The third exponent, rj s has proven to be numerically equal to S s = 5 thus 
ensuring that the hyperscaling law is fulfilled in a form appropriate for first order transitions and compact clusters 
p9|,p[. Moreover ,we have reported results of simulation for exponent S in case of NEKIM for the first time. 

Our present results give further evidence to the conclusion that the effect of fluctuations felt by the spin system at 
the PC transition is of interest in itself. 
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